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ABSTRACT 


An attempt is made to estimate aerodynamic parameters using the simulated flight 
data of various flight vehicles using Kalman filter technique. The present thesis 
demonstrates application of Kalman filter for the estimation of ballistic coefficient of a 
falling body by considering the effects of process noise. The mathematical model of a 
flexible aircraft involves larger number of stability derivatives as compared to a rigid 
aircraft. Furthermore, for a flexible aircraft, additional derivatives need to be included 
due to aeroelastic effects. Applicability of Kalman filter for parameter estimation is 
validated on simulated flight data generated for a rigid aircraft as well as for a flexible 
aircraft. It is concluded that Extended Kalman Filter method can be advantageously 
applied to estimate aerodynamic parameters firom flight data of a flexible aircraft. 
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CHAPTER 1 
Introduction 


All modern aerospace vehicles rely upon an understanding of dynamics and 
control to improve system performance. Successful system design requires an 
understanding of dynamic elements and the trade-off between vehicle dynamic 
characteristics, control system properties and system performance. Aircraft parameter 
estimation is one of the most outstanding and illustrated example of the system 
identification methodologies. The success of the system identification of the flight 
vehicle has been possible due to better measurement techniques and data processing 
capabilities provided by digital computers. Other factors that contribute to system 
identification are the developments in the fields such as estimation and control theory; the 
well understood basic principles of aerodynamic modeling and design of appropriate 
flight tests*’^’^. 

The equations of motion of early flight vehicles are derived fi'om Newtonian 
mechanics, usually assuming flight vehicles to be rigid. Out of the various forces and 
moments (aerodynamic, inertial, gravitational and propulsive) acting on a flight vehicle, 
the determination of aerodynamic forces and moments presents the most difficult 
challenge till date. Flight vehicle identification of high performance and highly 
augmented vehicles include nonlinearities affected by elastic modes, unsteady 
aerodynamics and erroneous data measurements. 

The concept of aircraft aerodynamic modeling was introduced by Bryan'^ , who 
provided relationship between the forces along the three Cartesian axes and the three 
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moments about these axes as a function of linear translational motion variables and 
rotational rates. The most commonly employed aerod 5 mamic models are typically 
formulated as truncated Taylor-series expansion of the form 

v-n dC: dC- 

■ i=X.Y,Z,l,m,n (1.1) 

where x„ are the aircraft motion variables (angle-of-attack, angle of sideslip, airspeed, 
angular rates, etc.) and U are the inputs (control surface deflection, etc.). X, Y, Z are the 
forces and 1, m, n are the moments about x, y and z axes respectively. However, such 
aerodynamic model structures can be surprisingly deceptive to their actual complexity. 
Some of the models postulated contain large number of parameters including two 
dimensional fimctions, acceleration derivatives, derivatives due to elastic modes, etc. 
The increase in the size of model structure makes the evaluation of the suitability of the 
model much more difficult. 

After a priori fix of the aerodynamic model, the next task of estimating 
parameters (stability and control derivatives) has been attempted by three different but 
complimentary techniques: analytical methods, wind tunnel methods, and flight - test 
methods. At the initial stages of aircraft design, analytical methods provide the only 
convenient way of estimating the aircraft parameters. However, the accuracy of such 
theoretical estimates being not sc high, there is need to verify these estimates with those 
obtained from wind tunnel testing and flight testing. In recent years, computational fluid 
dynamics (CFD) has contributed positively towards analytical estimation of aerodynamic 
coefficients by way of numerical solutions of total configuration via sophisticated and 
advanced Euler and Navier - Stokes flowsolvers^ . Nevertheless, experimental methods 
are essential to corroborate the analytical predictions. Although wind tunnel methods 


2 



improve the accuracy of estimation of parameters, they are time consuming and 
expensive. Furthermore, simulation of control surfaces, power effects and stringent flight 
conditions are difficult to simulate satisfactorily. Wind-tunnel estimates also suffer 
discrepancies due to interference effects of support systems, wall effects, turbulence 
level, etc. It is, therefore, desirable that the wind tunnel estimates be corroborated with 
the estimates from actual flight test data. 

One of the first attempts to obtain static and dynamic parameters from flight data 
was in 1947 by Milliken^ . Using semi- graphical method, he analyzed the frequency 
response data for obtaining the characteristics of the short period longitudinal motion of 
the aircraft. Later, Shinbrot^ proposed more general and rigorous methods based on 
application of ordinary and nonlinear least squares, for estimating parameters from 
transient maneuvers. In the late sixties and early seventies, the proliferation of high - 
speed, general purpose digital computers gave a new impetus to the development of more 
sophisticated, elegant and innovative identification techniques. The bulk of information 
on general identification theory and its dynamics have been documented in several 
reports and survey papers. A systematic discussion on the chronological evolution of 
flight vehicle system identification starting from determination of aircraft frequency and 
damping ratio from flight data in early twenties to the present day advanced applications 
is given in a survey paper of Hamel and Jategaonkar^ . 

The most commonly used methods to estimate stability and control derivatives from 
flight data are broadly classified into the following categories. 

• Equation error methods 

• Output error methods 
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• Filtering approach 

• Filter-error methods 

The principle of least squares* is used in the equation error method wherein the 
error gets minimized with respect to the unknown parameters in each of the equations. Its 
advantages include computational simplicity, non-iterative nature and applicability to 
both linear and non-linear models. The disadvantage is that the method cannot be directly 
applied if all the states are not measured accurately, and produces poor results if the 
measurements are noisy. Therefore in order to get accurate results, considerable effort is 
required for data reconstruction and smoothing. Sometimes these data processing tasks 
are more complicated than the task of parameter estimation itself. However these 
methods are useful as startup for more advanced estimators. 

In the output error method, the error between the measured and the model response 
is minimized. The method assumes that there exist no modeling errors. This method 
processes the measurement noise while assuming the model representation of the given 
system to be exact. A comprehensive survey of these methods is reported by Maine and 
Iliff^ . The methods like analog matching^ , Newton Raphson method*® , modified 
Newton Raphson method'® etc. fall in this category. Maximum likelihood (ML)^ method 
is one of the variants of output error methods. Maximum Likelihood estimates are those 
for which the observed value would be the most likely to occur; “most likely” is defined 
as maximization of likelihood function of the observed responses, given the parameters. 
The main advantage of this method is that parameter estimates are asymptotically 
unbiased, consistent and efficient provided the model assumed is correct. 



Kalman filter is a linear, optimal estimator of state variables of a linear, time- 
varying system, operating in a Gaussian stochastic environment. Filtering approach^’*^ is 
an extension of Kalman filter. Optimal estimator here is referred to a computational 
algorithm that processes measurements to deduce a minimum error covariance of the 
state of a system combining all the information available. 



Mathematical Model 

L ^ 

„ Augmentation of states 

by unknown parameters 

Estimation of augmented 
states by Extended 
Kalman Filter 

Fig. 1.1 Filtering Approach 

Generally, Kalman filter'^’*'^ is applicable and optimal for linear systems only. When 
either the system or the measurement equations are non-linear, the same algorithm can 
still be applied by local linearization of the system about the current state. Such filter 
applied to nonlinear systems is called Extended Kalman Filter (EKF) , and it need not be 
optimal. EKF produces estimates of the parameters that approximately mmimize the 
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mean square error in the parameter estimates themselves as opposed to ML and least 
squares which minimize a cost function that is based on matching the output variable 
behavior given a specific input trajectory. The filtering approach as presented pictorially 
in Fig. 1 . 1 is well suited for online applications and can handle process noise (modeling 
errors) and measurement noise (instrument errors). Extended Kalman Filter produces 
time histories for each of the estimated parameters, which is very useful if the time- 
varying nature of the parameters is of interest. But the performance of the filter strongly 
depends on the initial statistics of process noise, measurement noise and a priori 
covariance matrix. Information about the measurement noise covariance can be obtained 
firom the laboratory calibration of measurement sensors used. Process noise is however 
more difficult to determine. Several procedures for obtaining measurement and process 
noise prior to estimation are available in literature. Well known methods by Mehra''* and 
Morelli*^ are being used with partial success. Other kind of approach known as adaptive 
filtering adjusts the variances during the estimation of parameters. Myers and Tapley*^ 
gave one such algorithm of adaptive filtering/tuning. In the present work no attempt has 
been made to use adaptive filtering technique. 

Filter error method^’*^ is the most general stochastic approach to parameter 
estimation, which accounts for both process and measurement noise. It can give good 
estimates from data collected in turbulent atmosphere (process noise) also. In this direct 
approach an appropriately defined cost function, usually the statistically formulated 
likelihood function, is optimized with respect to the unknown parameters using some 
suitable optimization method. The Gauss Newton method is generally used for this 
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purpose. Filter is used only for the natural purpose of obtaining the true state variables 
from the noisy measurements, i.e. for state estimation. 

Parameter estimation of an aircraft with many measurements like, acceleration 
(both linear and angular), angular orientation, speed, angle of attack etc. has been 
discussed in detail in Ref,. 5. But from cost effectiveness point of view, it may not be 
feasible to use many sensors for air borne vehicles that go through man y development 
trials. 

In the present work, EKF method is used for parameter estimation from flight data 
of air borne vehicles. The motivation here was to conduct a study on the applicability of 
the method in extracting aerodynamic parameters by processing flight data obtained for 
different class of flight vehicle/store. The method has been applied to flight data of a one- 
dimensional ballistic target, flight data of a rigid aircraft in longitudinal short period 
mode and also to flight data of a flexible aircraft in longitudinal short period mode. Any 
parameter estimation method requires adequate information about vehicle dynamics to 
estimate aerodynamic parameters correctly. Flight data of one-dimensional ballistic target 
contains very limited information regarding its motion/control variables. The flight data 
obtained through aircraft maneuvers contains more information regarding its motion and 
control variables. The applicability of the EBCF is tested for these three different classes of 
fli gh t data. It is observed that EKF can be advantageously applied on flight data of 
flexible aircraft to estimate few aerodynamic parameters with acceptable level of 
accuracy. 

Chapter 2 describes the description of Kalman filter and the mathematics involved 
in the filter design. Chapter 3 presents methodologies used to generate flight data for 
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different type of vehicles considered. In chapter 4 results and discussions are presented 
and chapter 5 concludes with scope for future work. 
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CHAPTER 2 


Kalman Filter 


Kalman filter is an optimal recursive data processing algorithm and it 
incorporates all information that is available to the filter. It processes all available 
measurements, regardless of their precision, to estimate the current values of the 
variables of interest with use of 

knowledge of the system and measurement device dynamics 

the statistical description of the system noises, measurement errors, and 

uncertainty in the dynamic models 

- any available information about initial conditions of the variables of interest 

This filter performs the conditional probability density propagation for 
problems in which the system can be described through a linear model and in which 
system and measurement noises are white and Gaussian. The three basic assumptions 
in Kalman filter formulation are 

- the model is considered to be linear . 

- the measurement noise and system noise are white. 

- the measurement noise and system noise are Gaussian. 

The physical implications of these assumptions are discussed. A linear model 
is justifiable because when nonlinearities do exist, the typical engineering approach 
is to linearize about some nominal point or trajectory, achieving a perturbation 
model or error model. 

“Whiteness” implies that the noise value is not correlated in time. The 
knowledge of the present noise value is no way helpful to predict the noise value of 
any other time. Whiteness also implies that the noise has equal power at all 
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frequencies. A system will be driven by a wideband noise, noise having power at 
all frequencies above the system bandpass and essentially constant power at 
frequencies within the system bandpass. The assumption of noise as white would 
extend this constant power level out across all frequencies. 



Fig. 2.1 Power spectral density bandwidths 
Within the bandpass of the system of interest, the fictitious white noise looks 
identical to the real wideband noise. The mathematics involved in the filter are vastly 
simplified by replacing the real wideband noise with a white noise which from the 
system’s “point of view” is identical. Gaussianness is related with amplitude whereas 
whiteness is related with frequency. The probability density of Gaussian noise 
amplitude takes on the shape of a normal bell -shaped curve. This assumption is 
justified physically by the fact that a system or measurement noise is typically caused 
by a number of small sources. Mathematically, when a number of independent 
random variables are added together, the summed effect can be described very closely 
by a Gaussian probability density, regardless of the shape of the individual densities. 
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The first and second order statistics (mean and variance or standard deviation) 
of a noise process can be easily known. In the absence of any higher order statistics, 
there is no better form to assume than the Gaussian density. The first and second order 
statistics completely determine a Gaussian density, unlike most densities which 
require endless number of orders of statistics to specify their shape entirely. The 
Kalman filter, which propagates the first and second order statistics, includes all 
information contained in the conditional probability density. 

The real world situation is described by a set of nonlinear differential 
equations to apply extended Kalman filter techniques. These equations are expressed 
in nonlinear state-space form as a set of first-order non linear differential equations as 

x = f{x)+w (2.1) 

where ‘ x’ represents system space, ‘/(^)’ is a nonlinear function of those states 
and ‘w’ is a random zero-mean process. The continuous process-noise matrix 
describing the random process w for the preceding model is given by 

[j2] = jE(wM'^) (2.2) 

The measurement equation, required for the application of extended Kalman filtering, 
is considered to be a nonlinear function of the states according to the equation 

z =h{x)+v (2.3) 

where ‘v ’ is a zero-mean random process described by the measurement noise matrix 
[R], which is defined as 

[ii]=£(w’-) (2.4) 

For systems in which the measurements are discrete, the nonlinear measurement 
equation is written as 

z*=/i(x*)+v* (2.5) 
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The discrete measurement noise matrix [i?]k consists of a matrix of variances 
representing each measurement noise source. Since, the system equation, Eq. 2.1 and 
measurement equation, Eq 2.3 are nonlinear, a first order approximation is used in the 


continuous Riccati equations for the manipulation of systems dynamics matrix [f] 
and the measurement matrix [if] . The matrices are related to the nonlinear system 
and measurement equations according to the relations 





( 2 . 6 ) 




(2.7) 


The fundamental matrix [^^], required for the discrete Ricatti equations, can 
be approximated by the Taylor-series expansion for exp([F] Ts) and is given by the 
equation 


[4=W+M3; 



3! 


( 2 . 8 ) 


where ‘Ts’ is the sampling time and [/] is the identity matrix. In our applications of 
extended Kalman filtering, the series is approximated by only the first two terms, 
because \_4'\k only used for the calculation of Kalman gains and the matrix may 
not necessarily improve the performance of the filter by considering more terms. 
Therefore [(f>\ is given by 

W (2-«) 
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The matrix Ricatti equations, required for the computation of the Kalman 
gains, are given by the equations 



( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


where is the covariance matrix representing errors in the state estimates after an 
update, [iiTi] is the Kalman gain and [Af] is the covariance matrix representing errors 
in the state estimates before an update. The discrete process-noise matrix \q]^ can be 
found from the continuous process-noise matrix according to equation 

'G],‘]W)MW)Y dr ( 2 . 13 ) 

0 

The preceding approximations for the fundamental and measurement 
matrices are used in the computation of the Kalman gains. The new state estimate 

is the old state estimate projected forward to the new sampling x,^ plus a gain 
times a residual. The residual is the difference between the actual measurement 
and the nonlinear measurement ) . 

x,^x,+Kdz,-h{x,)] (2.14) 

The old estimates that have to be propagated forward do not have to be done 
with the fundamental matrix but instead can be propagated directly integrating the 
actual nonlinear differential equations forward at each sampling interval. Euler 
integration is applied to the nonlinear system of differential equations and is given by 
the equation 

Xt=x,,_i+x^_J, (2.15) 
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where the derivative is obtained from 


V,=/(V.) ( 2 - 16 ) 

In the preceding equation the sampling time Ts is used as an integration 
interval. In the problems where the sampling time is large, Ts would have to be 
replaced by a small integration interval, or possibly a more accurate method of 
integration has to be used. 

Derivation of Scalar Riccati Equations 

When there is no deterministic disturbance or control scalar the discrete model 
is given by 

Xt=<l>kXk-i+^k (2-17) 

where ‘ ’ is the scalar that propagates the states from one sampling instant to the 

next and ‘ w* ’ is the white process noise. Hie scalar Kalman - filtering equation for 
this discrete model is 

x,=4>kXk-^^K,{zk-Hc|>kX,.,) (2.18) 

The measurement scalar equation is given by 

z,=Hx,+v, (2.19) 

where ‘ ’ is the measurement scalar that relates the state to measurement and ‘ v* ’ is 

the measurement noise. The error in the estimate at time ‘A’ is given by this equation 

Xk =Xt -x^ =x^-<j>kXk-x +K„ {z^ ( 2 . 20 ) 

Substituting Eq. 2.17 and Eq. 2.19 in the Eq. 2.20 yields the following equation 

Xk Xk-i +'^k -AXk-i -Kk W^x^.i +Hwt +V^ -H^kXk-i ) (2-21) 
■ The error at time ‘ k-1 ‘ is given by the equation 
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( 2 . 22 ) 


Combining similar terms of the Eq. 2.21 yields 


x,={\- +(1 - -K, 


(2.23) 


Pk , Qk and Rk are defined as covariance scalar, process noise scalar and 
measurement noise scalar respectively and are computed by the expectations of 
square of the error in the estimate, square of the process noise and square of the 


measurement noise respectively as follows 


Pk=E{x,) 


(2.24) 


Qk=Ei^l) 


(2.25) 


Rk=Eivl) 


(2.26) 


Squaring and taking expectations of the Eq. 2.23 yields 

P, = (1 - {P,Jl +Q,) + Kl R, 


( 221 ) 


since the process noise and measurement noise are considered to be zero mean noise 
their expectations are zero. 

The Eq. 2.27 can be further simplified by defining a term Afias 


M,=P,Jl+Qk 


(2.28) 


Eq. 2.28 is analogous to the Ricatti first equation, Eq. 2.10 and on substitution of Eq. 


2.28 in Eq. 2.27 the equation simplifies to 


P,={\-K,HfM,+KlR, 


(2.29) 


The optimal gain that will minimize the variance of the error in the estimate is 
obtained by equating the derivative of Pk with respect to Kalman gain, Kk to zero. 

.^=0 ( 2 . 30 ) 


2(1 - K^H)Mi,i-H)+2K,R, =0 


(2.31) 
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The Kalman gain, is given by the expression 




(2.32) 


Eq. 2.32 is analogous to the second Riccati equation, Eq. 2.11. Substitution of the 
optimal gain Kk in the Eq. 2.29 yields. 




1 — 


^ ^ MM 




Mk + 


H^M,+R, 


R. 


(2.33) 


The Eq. 2.33 is simplified further by the use of Eq. 2.32 to 

p = _ ^kKk ( 2 . 34 ) 

* H 

Rearranging the optimal gain equation, Eq. 2.32 

KtRk=MM-H''MMk (2.35) 

Substituting Eq. 2.35 in Eq. 2.34 yields 

( 2 . 36 ) 

H 

Eq. 2.36 is analogous to the third Riccati equation, Eq. 2.12.These three Riccati 
equations are simply an iterative way of finding the optimal gain at each time step. 
The matrix Riccati equations also resemble the scalar Riccati equations and their 
derivation is far more complex than the scalar derivation and is not mentioned in the 
present work. 
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CHAPTER 3 


Generation of Flight Data 

Due to the non-availability of real flight data, simulated flight data is 
generated for different class of flight vehicles for the purpose of parameter estimation. 
Mathematical models used to generate flight data of ballistic target (FD-BT), longitudinal 
motion of rigid aircraft (FD-RA/C) and longitudinal motion of flexible aircraft (FD-AE 
A/Cl, FD-AE A/C2) are presented in subsequent paragraphs. For the generation of 
simulated flight data of ballistic target the acceleration equation governing the path of the 
trajectory was solved using Euler integration. In the case of aircraft, longitudinal 
equations of motion were solved using fourth-order Runge -Kutta method. Subsequently, 
the simulated flight-data of the air borne vehicles is used in Kalman filter algorithm and 
processed for the estimation of aerodynamic parameters. The equations of motion are 
needed for generating the simulated flight data and also in the Kalman filter algorithm for 
the propagation of states from one time step to the other. 

Flight Data of Ballistic Target (FD-BT)* 

Knowledge of target ballistic coefficient is used in advance guidance laws such as 
predictive guidance to relax the interceptor acceleration requirements. In addition, 
knowledge of the target ballistic coefficient is required for fire control due to the 
importance of accurate intercept point predictions in launching the interceptor on a 
collision course. Therefore, accurate estimation of ballistic coefficient of a target re- 
entering the atmosphere is very important for both guidance and fire control purposes. 
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The flight data for simulating such a vehicle motion is modeled to investigate the 
applicability of EKF method^ in extracting parameter (ballistic coefficient) from flight 
data. 

The one-dimensional example of a ballistic target failing on tracking radar is 
considered. The target was initially at 2, 00,000 ft above the radar and had a velocity of 
6000 ft/s towards the radar, which is located in the surface of a flat Earth. The trajectory 
of the ballistic target is presented in Fig. 3.2. The radar measures the altitude of the target 
with 25-ft standard deviation measurement accuracy. The radar picks measurements for 
every 0.1 -sec. The simulation is done for 30 sec. An extended Kalman filter is built to 
estimate the altitude, velocity and ballistic coefficient. 


Drag 

(Unknown) 


6000 Ft/Sec 


9 

(Known) 


200,000 Ft 


Radar 

Fig. 3.1 Radar tracking of ballistic target . 

The two forces acting on the object are the drag force and the gravity force. The 
equation that governs the motion of the object is given by 
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(3.1) 



where ‘ x ’ is the acceleration acting on the object gg ’ is the acceleration due to gravity 
and ‘j2p’ the dynamic pressure and is given by the equation 

Q,=0.5pgx^ _ (3.2) 


where ‘x’ is the velocity of the target and the air density ‘Pg ’ in Eq. (3.2) is an 
exponential function of altitude and is given by the equation 

Pg =0.00346'^^^^““ (3.3) 

The term ‘ yff ’ in Eq. (3.1) is the ballistic coefficient and is expressed as 


W 

13= ® 


^ref^DO 


(3.4) 


where ‘JTb ’ is the weight of the ballistic target, ‘ ‘ is the reference area of the object 

and ‘ Cqo ’ is the coefficient of drag. The acceleration acting on the object is expressed as 
the nonlinear second order differential equation as 


.. _^034gge 
X 


-JC/ 22000 


2y9 


■ga 


(3.5) 
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Fig. 3.2 Ballistic target trajectory for different values 
The procedure of application of Kalman filter algorithm for the case of Ballistic Target is 
explained in Appendix -A. 

Simulated Flight Data of Rigid Aircraft(FD-RA/C): 

In most cases, longitudinal maneuvers predominantly excite tiie short period 
mode and not the phugoid mode. In short period mode the flight velocity is essentially 
constant during the maneuver. 3-2-Tl control (elevator) input is considered for the 
excitation of short period mode. The example aircraft mentioned in Ref. 17 was chosen 
for the present study. The mass, geometric and aerodynamic characteristics of this 
example aircraft are presented in Table 3.1. Trim flight conditions correspond to straight 
and level cruise flight at an altitude of 1500m and at a Mach number of 0.6. The rigid 
body short period longitudinal response to a given elevator input was simulated for 8 
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seconds. The short period longitudinal response was simulated using the following 
equations. 

d-q = - puSCj^Ilm (3.6) 

q = pu^ScCj2I, (3.7) 

The equations that define coefficient of lift {Cl ) and pitching moment coefficient 
(Cm ) in Eq. (3.6) and Eq. (3.7) to describe the aerodynamic model are presented in Eq. 
(3.8) and Eq. (3.9) 

Cl — Clq + C^^GT + — + Cls^ Sg (3 .8) 

C.=C„+C.,a + C„^+C.,S, (3.9) 

The pitching moment coefficient {€„) is vfith reference to the center of gravity . 
The Eq. (3.6) and Eq. (3.7) were solved using fourth-order Rimge-Kutta algorithm with a 
time step of 0.001 sec to obtain simulated fli^t data. The contribution due to aeroelastic 
effect was neglected for rigid aircraft case. The flight data used for parameter estimation 
is pictorially presented in Fig. 3.3 and Fig. 3.4. This flight data will be referred to as FD- 
RA/C for further purpose. The elevator input (3-2-1-1) used to generate this flight data is 
shown in Fig. 3.5. 


21 


. 1 . 5219 - 2 — 








Table 3.1 Geometric, mass and aerodynamic characteristics of example aircraft 


Mean Chord 

c = 4.663 m 

Wing Sweep 

A=65° 

Wing Area 

180.79 m^ 

Mass 

m= 130642.31 kg 

Moment of Inertia about Y axis 

/y= 8677231.17 kg-m^ 

Ml = 248.86 kg-m^ 

Modal generalized masses 

M 2 = 12984.23 kg-m^ 

Ms = 1807.34 kgW 

M4 = 59048.57 kgW 

Aerodynamic Parameters 

Clo 

0.34 

Cla 

2.9223 

Clq 

-14.7 

n 

0.43548 


CjnO 

-0.252 

Cma 

-1.6617 

^tnq 

-34.75 


-2.578 
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Table 3.2 Vehicle Configuration 


In vacuo modal frequencies <y, (rad/s) 


Configuration 

Mode 1 

Mode 2 

Mode 3 

Mode4 

AEA/Cl- 

12.57 

14.07 

21.17 

22.05 

AEA/C2 

6.29 

7.04 

10.59 

11.03 


Simulated Flight Data of Flexible Aircraft FD-AE A/Cl & FD-AE 


Aid) 

The non linear equations of motion for a flexible aircraft are considered which are 
obtained from Ref. 17. 


d - q - puS jlm 


Cio + + Ci^g^ 5^ 


+ 




i=l 


(3.10) 


q = pu 


■Sc Illy 


^mq 


qcjlu + C„sS^ + ^ (c„^. 77 ,. + 7 ,c/ 2 w) 


(3.11) 


L /=! J 

v/here angle of attack a, pitch rate q and control input 5^ represent the 
perturbations and 77 , and 7 , are the generalized displacement coordinates and their 
derivatives. ,C,^, Ci^g^, C„„, C„,, C^g^, and C„^, are the 

aerodynamic coefficients as defined in Ref. 17. The second order differential equation 
that is satisfied by the generalized displacement coordinates and its derivatives is taken 
from Ref. 17 and is given in Eq. 3.12. The additional term 2 ^, 67 , 77 ^ representing the 

structural damping is also included. 
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Vi +^^i^A +^fVi = pu^ScIlM, 


c*a+c;' qci2u*cfs*f}pin, +q; njcjiu) 

J=i 


(3.12) 


where 6?, , and Mi are the in-vacuo frequency, modal damp ings and modal 


generalized mass respectively and CJ' , C^' , CJ' , and C|;' represent the generalized 

force derivatives due to coupling in elastic and aerodynamic degrees of freedom. 
Eq. 3.10, Eq. 3.11 and Eq. 3.12 are integrated to generate motion variables a and^ . The 
elevator control input is 3-2-1-1 and is presented in Fig. 3.5. The aerodynamic model 
incorporating aeroelastic effects is presented in Eq.(3.13) and Eq.(3.14) 


/=! 

f +(Q,,7,)c/2m] 

(3.13) 

l=:l 

li +[c„f„v)cl2u] 

(3.14) 


The nmnerical values of the aerodynamic parameters are listed in Table 3 .3. The 
flight data obtained for this case is pictorially presented in Fig. 3.6 and Fig. 3.7. To study 
the effect of flexibility of the structure on parameter estimation, two cases are considered. 
The two set of modal frequencies for the two different configurations are listed in Table 
3.2. To start with aircraft with moderate flexibility is considered. The flight data obtained 
with this case will be referred to as FD-AE A/C. In the second case more flexibility is 
included and the flight data generated for this case is referred to as FD-AE A/C2. 
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Table 3.3 Aerodynamic Parameters of the flexible aircraft from Ref. 17 





tTi 

VO 



VO 

CO 

VO 

vn 

VO 

j 

1 

CM i 


<u 


o 

o 

o 

O 

O 

o 

O 

o 

o 

O 1 

o 

VO 

O 

<D 

0) 

1 

0) 

o 

1 

o 

1 

a> 

<!> 

1 

o 

<i> 

<l> 

d 


d 

m 

VO 


04 

00 

0^ 

trv 

Os 

VO 

ON 


m 


oo 

tn 

04 

CO 

VO 


vn 

O 

VO 

ON 


cn 


<N 

o4 

oi 

1 

OS 


vn 

i 



K 




04 



CO 




04 

VO 

CO 


<D 


o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

ON 

o 

<D 

1 

<u 

1 

<D 

1 

a> 

i 

i 

4> 

a> 

t 

<u 

<D 

i 

d 

oo 

o 

0- 

VO 



VO 


On 

»o 

tO' 

r— ( 




Os 


Os 


rj- 


CM 

ON 

t> 

CO 




cn 


04 

1— H 

CO 

1 

»— 1 

VO 

1 


OO 

1 

1 


CO 


CO 

CO 1 

o 

o 

1 

o 

o ■ 

i> 

<u 

(D 

<D 

r—< 

oo 

On 

to 


On 


T— ( 

00 

CM 

1 




o o 




tn 00 
CO 

m CN 


CN 


O 9 


Os CS ^ 

fs ^ ^ 

9 o o 


E E EE BE I E I 

oooooooou 


► s 6 
u o 


(S 

(N O 


;C| g S 2 § “ 

2 2 8 3 o S " 

I o o I ? o d 


00 Os 

O o 
o Q 


•S’ 2 5 *5 

o o u o 


5 3 

O o o 


cs 


o O I o 
d d I o 


o o 
d d 


o O 


S? :r ^ 

o o o o o 


s? ^ 

u o 


18 





































































CHAPTER 4 


Results and Discussion 

In this chapter, estimated parameters obtained through EKF method are presented. 
In the case of ballistic target the factors that affect the behavior of the filter are discussed. 
The simulated data corresponding to ballistic target (FD-BT), rigid aircraft (FD-RA) and 
flexible aircraft (FD-AE ACl & FD-AE AC2) are used as the measured flight data. The 
estimated parameters are presented along with their standard deviations to assess the 
accuracy of the estimates of aerodynamic parameters. 

Estimation of Ballistic Coefficient from Flight Data of Ballistic Target 
(FD-BT) 

As mentioned in Chapter 3 the ballistic target is assumed to fall on a straight line 
path towards the surface based tracking radar. The balhstic coefficient (y0 ) is estimated 
through EKF method. It has been observed that there is a negligible change in the 
estimated value of p by the increase of terms in the Taylor series expansion for the 

approximation of fundamental matrix [a] in Eq. A. 10 of Appendix A. This is valid 
because the fundamental matrix is actually an infinit e Taylor series expansion of the 
product of sampling time ts and system dynamics matrix [f]. 

Table 4.1 Estimated ballistic coefficient without process noise 



True Value 

Estimated Value by 

2 terms 

Estimated Value by 

3 terms 

Ballistic Coefficient 
p (Lb/Ft^sec^) 

L-_- — 

500 

497.8827 

(0.290904)* 

499.8196 

(0.290046)* 


Standard Deviation 
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For the same case the process noise is included assuming that the filter’s 
knowledge of the real world is in error. It has been observed that addition of process 
noise to the filter increases the errors in the estimates. Also there is not much difference 
by increase of number of terms in the Taylor series expansion for the approximation of 
fundamental matrix From the Table. 4.1 and Table. 4.2 it is evident that by 

including process noise the accuracy in the estimates deteriorated. Similar observation 
can be made by comparing Fig. 4.1 and Fig. 4.2 with Fig. 4.3 and Fig. 4.4. 

Table 4,2 Estimation of ballistic coefficient with process noise 



True Value 

Estimated Value by 

2 terms 

Estimated Value by 

3 terms 

Ballistic Coefficient 
fi ( Lb/Ft/sec^ ) 

500 

488.9278 

(11.07214) 

488.7636 

(11.23639) 



Time . sec 


Fig. 4.3 Estimation of ballistic coefficient ji with process noise 
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Fig. 4.4 Error in the estimation of ballistic coefficient /? with process noise 


Estimation of Aerodynamic Parameters from Flight Data of Rigid 
Aircraft (FD-RA/C) 

Knowledge of aerodynamic parameters is of paramoimt importance, to develop 
accurate mathematical model which represents the longitudinal dynamics of an aircraft. 
The aerodynamic parameters needed to develop the aerodynamic model are force 
derivatives and pitching moment derivatives Thus, a 

study was carried out to explore the possibility of extracting the parameters using EKF 
method. The simulated data generated, FD-RA/C of the example aircraft is added with 
Gaussian noise and is used as the measured data in the EKF algorithm. The six 
aerodynamic parametersQ„ and C^^ are considered as state 

variables along with the flight variables, a and The algorithm is run for an elevator 
input 3-2-1-1 for 8 seconds and the estimated aerodynamic parameters are tabulated 
along with their standard deviations in Table. 4.3. The measurement noise variance is 
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0.016 for both a and q. In Case A the values of diagonal elements of process noise 
matrix are same and are randomly chosen to be 0.008 and in Case B the values are 
differently chosen as an attempt to tune process noise. 

Table. 4.3 Estimated Aerodynamic Parameters from FD-RA/C 



Aerodynamic Parameters 

Case 

Q. 


c 

'^LS, 

Cma 

^mq 


True Values 

2.922 

-14.7 

0.435 

-1.66 

-34.75 

-2.578 

No Process Noise 

2.9096 

(0.00764) 

-22.7217 

(0.75708) 

0.6756 

(0.01976) 

-1.6078 

(0.00166) 

-38.3583 

(6.68809) 

-2.5619 

(0.00289) 

Process Noise 

Case A 

2.8772 

(0.06228) 

-16.6971 

(1.78776) 

-0.3212 

(0.13688) 

-1.8501 

(0.03303) 

-39.0636 

(6.9463) 

-2.7122 

(0.10493) 

Process Noise 

Case B 

2.8887 

(0.08767) 

-16.8602 

(1.89439) 

-0.2328 

(0.14417) 

-1.6099 

(0.0784) 

-41.527 

(6.94662) 

-2.8064 

(0.0588) 






Fig . 4.6 Comparison of the estimated q value with that of FD-RA/C. 


Through the values listed in Table 4.3 it can be observed that by including process 
noise the estimated parameters show large variance however we can also infer that by 
proper choice of noise values the estimate of parameters can be improved. 

Estimation of Equivalent Aerodynamic Parameters from Flight Data of 

Flexible Aircraft (FD-AE A/Cl & FD-AE A/C2) 

The simulated data from FD-AE A/Cl and FD-AE A/C2 is used in the EKF 
algorithm by adding Gaussian noise to the data. All the 4 elastic modes of the test aircraft 
are considered in the generation of the simulated data. A full order model of an 
aeroelastic aircraft has too many parameters to yield satisfactory estimates using any of 
the conventional parameter estimation methods. In view of this, a study was earned out to 
identify a simplified model with reduced number of parameters, and to evaluate how the 
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resulting parameters are affected by model simplifications. To start with, a rigid body 
model was assumed and parameters were estimated from flight data that contain the 
aeroelastic effects. It is expected that the parameter estimates thus obtained would absorb 
the aeroelastic effects. For the convenience of discussion, these parameters are referred 
by equivalent parameters’ in Ref 19. An analytical expression has been proposed to 
anal 3 dically compute the numerical values of the equivalent parameters. It is shown that 
the numerical value of analytical expression indicates the degree of flexibility of the 
aircraft, and thereby, a criterion based on it is suggested for deciding adequacy or 
otherwise of using simpler rigid body models in estimation algorithm. In the EKF 
algoritlim only aerodynamic parametersQ„,C^,Q^^,C„^,C„^ and are estimated 
and presented in Table.. 4.4. 


Table 4.4 Equivalent Aerodynamic Parameters from FD-AE A/Cl and FD-AE A/C2 



Equivalent Aerodynamic Parameters 

Q. 



^ma 

^mq 


Tme Values 

2.922 

-14.7 

0.435 

-1.66 

-34.75 

-2.578 

FD-AE A/Cl 

2.4927 

-1.5588 

0.9644 

-1.3141 

-20.6907 

-2.2514 


(0.006445) 

(0.6973) 

(0.0179) 

(0.0012) 

(6.6134) 

(0.00271) 

FD-AE A/C2 

1.7955 

-7.7374 

-0.2195 

-0.4582 

-15.2807 

-1.5197 


(0.017004) 

(2.56618) 

(0.06234) 

(0.001378) 

(6.932) 

(0.00752) 


The equivalent aerodynamic parameters are in close agreement with values presented in 
Ref 19. For the case of FD-AE A/C 1 the measurement noise for a and ^ is 0.016 but 
for the case of FD-AE A/C 2 the measurement noise is 0.09. For the case of FD-AE A/C 
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2 the aeroelastic effects of more flexible aircraft are got included in the algorithm as 
measurement noise. 

Estimation of Aerodynamic Parameters from Flight Data of Flexible 
Aircraft (FD-AE A/Cl & FD-AE A/C2) 

In the final case both the flight data FD-AE A/Cl and FD-AE A/C2 are processed 
by the EiCF algorithm. The algorithm includes only 15 parameters, taking into account 
only the first elastic mode. Along with the flight variables a and q the generalized 
displacement coordinate and its derivative are also taken as the state variables. The 
estimated values of aerodynamic parameters are presented in Table. 4.5. It is seen that the 
estimation deteriorates as the flexibility increases. Thus for high flexibility aircraft such 
an approximation may not yield better results Comparison of estimated response of a ,q 
with that of the flight data FD-AE ACl and FD-AE AC2 are presented pictorially in 
Fig.4.7 to Fig. 4.10. Based on these Fig. 4.7 to Fig. 4.10 it can be easily seen that the 
estimated response {cc,q) closely matches with the simulated response. However, there is 
a large difference between these parameters; it can be used for simulators and control law 
specifications to initiate the analysis. 
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Table 4.5 Estimated Aerodynamic Parameters from flight data FD-AE ACl and 
FD-AE AC2 



True Values 

FD-AE ACl 

FD-AE AC2 


2.922 

3.2861 

(0.0445) 

3.0798 

(0.03598) 

O 

-14.7 

-13.2036 

(1.02) 

-24.6586 

(1.7327) 

Oq 

0.435 


0.2619 

(0.03821) 


-1.66 




-34.75 

-32.4563 

(0.2808) 

-41.1266 

(0.49461) 


-2.578 

-2.565 

(0.0145) 

-2.1175 

(0.01376) 

Q 

0.0288 

0.0556 

(0.0024) 

0.012 

(0.000479) 

Q 

h\ 




'”171 

0.0025 

-0.0411 

(0.0025) 

-0.0201 

(0.000129) 

C.., ^ 

-0.159 

-0.0799 

(0.022) 


C,.. 

-0.014898 

-0.017 

(0.0018) 

-0.0244 

(0.000281) 

Ql 

Vi^q 

-0.0949 

-0.1 

(0.0186) 

-0.0672 

(0.009936) 

q; 




q! 

-0.00042 


lIHliHH 

q; 

-0.012835 

■llinQQiHi 

HHH 
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CHAPTER 5 


Conclusion and Suggestions For Future Work 

5.1 Conclusion 

In the present work, EKF method has been applied to estimate aerodynamic parameters 
from simulated flight data. The method has been applied starting from flight data of a 
one- dimensional ballistic target to flight data of a rigid aircraft and then to the flight data 
of a flexible aircraft. It is observed that EBCF method can be applied successfully to 
estimate the parameters from the flight data of the three cases. The method also estimates 
equivalent aerodynamic parameters in which aeroelastic effects get absorbed. If 
appropriate mathematical model of the system is provided then the method can be used 
advantageously to estimate force and moment derivatives. 

5.2 Suggestions for future work 

1. Due to non-availability of real flight data, aerodynamic parameters were estimated 
using simulated flight data. In order to enhance the proposed schemes, these parameters 
should be validated using real flight data. 

2. All the simulated data is generated under standard atmospheric conditions. 
Applicability of the proposed schemes to estimate parameters under non-standard 
atmospheric conditions can be explored. 

3. To obtain reliable results from recursive techniques like EKF, it becomes necessary to 
tune the filter parameters fairly well offline before implementing the techmque online. 
Proper filter tuning methods like Tapley Meyer’s algorithm^^ or Adaptive Kalman 
filtering*^ can be applied to tune the process noise and measurement noise characteristics. 
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APPENDIX -A 


Kalman filter design for Ballistic Target 

The ballistic coefficient of the object is xmknown and it is made one of the states and also 
it is assumed that the ballistic coefficient is constant. The states for the proposed filter 
are given by 


[X]= 


(A.1) 


The complete set of equations required to broild the filter are 


X = 


0.0034gj5e 


-x/ 22000 


IP 


'Sb 


(A.2) 


y0=O (A.3) 

The resultant filter can be made more robust by adding process noise to the derivative of 
the ballistic coefficient. Then the derivative of the ballistic coefficient becomes 

,8-u, (A.4) 

The two differential equations Eq. (A.2) and Eq. (A.3) can be linearised as follows 


Ax 

Ax 

Ay5 


dx dx dx 
dx dx dp 
fix Sx Sc 
dx dx dp 

M M M 

dx dx dp 


Ax 


’0 ■ 

Ax 

+ 

0 





(A.5) 


For the present case it is assumed that the ballistic coefficient is constant and so Wj is 
zero. 
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The systems dynamics matrix is the matrix of partial derivatives and is given by 


[F] 


dx 

dx 

dx 

dx 

dx 


dx 

dx 

d& 

dx 

dx 


M 

M 

d'p 

dx 

dx 

dp 


(A.6) 


The partial derivatives are evaluated at the current state estimates. The system dynamics 
matrix of the present case turns out to be 


[F] = 


0 1 0 

-PbSb^^ PbSb^ -PbSb^^ 

44,000y? 

0 0 0 


(A.7) 


For simplification the following terms are introduced 


r _ -PbSb^ 
J2\ ~ 


44,OOOy0 
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_ PbSb^ 


fn ~ 


-PbSb^ 

10^ 


The system dynamics matrix transforms to 


[F] = 


0 1 0 

fi\ fn fn 
0 0 0 


(A.8.1) 


(A.8.2) 


(A.8.3) 


(A.9) 


The fundamental matrix is found by the Taylor-series approximation and is 
given by the equation 
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m)] = [/] +[F]t + 


+ . . . 


(A.10) 


Mll+M?! 

2! 3! 


Since the fundamental matrix is only used in the Ricatti equations the first two terms of 
the Taylor-series expansion is considered to obtain the fundamental matrix. Also by the 
increase in niimber of terms there may not necessarily be a considerable change in the 
estimate of the states. 


[(^(0] « m HF]t 


1 t 0 

0 0 1 


(A.11) 


The discrete fundamental matrix is obtained by substituting Ts for t and is given by 


m 


1 T, 0 

flX^s fll^s 

0 0 1 


(A. 12) 


In the example considered the altitude measurement is a linear function of the states and 
is given by 



(A.13) 


where v is the measurement noise. By inspection of the measurement equation we can 
get the measurement matrix [i/] and the discrete measurement noise matrix [R] 


[H] = [l 0 O] (A. 14) 

R, = £(v,v[) (A.15) 


The discrete measurement noise matrix for the problem considered turns out to be a 
scalar and is given by the equation 

<A.16) 


A2> 



If the resultant filter needs to be more robust then process noise can be included. The 
continuous process noise matrix [Q] is given by 


[Q] = 


0 0 0 
0 0 0 
0 0 <!,, 


(A.17) 


where is the spectral density of the white noise source assumed to be on the 

derivative of the ballistic coefficient .The discrete process noise matrix [Q\ is derived 
from the continuous process-noise matrix according to the expression 


Kj = (A. 18 ) 

0 

After integration the final expression for the discrete process noise matrix is given by the 
expression 


[a] = 


0 


0 


t! 


/b 

/23^ 


3 

2 


f. 


23 


The Kalman-filtering equations are written as 


(p^l ~ ^k) 

4 = ^ + ^ 2 * (4 -^k) 

Pk ~ Pk-\ •^3t i^k ~ ^k) 


(A.19) 


(A.20.1) 

(A.20.2) 

(A.20.3) 


The three-state extended Kalman filter and Riccati equations were programmed, along 
with the simulation of the real world. The process noise matrix is set to zero assuming 
to be a constant. The third element of the initial covariance matrix is set to 300^ some 
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guess value to reflect the uncertainty in the initial guess at the estimated ballistic 
coefficient. The initial states of the filter are chosen close to the true values. The altitude 
state X is in an error by 25 ft and the velocity state x is in an error by 150 ft/s. The initial 
covariance matrix reflects the errors in the position state and the velocity state. 
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